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1. Introduction 

The study of hadrons allows experimental access to a wide variety of phenomena. However, at 
low energies QCD does not provide easy theoretical tools for parsing experimental results into mea- 
surements of fundamental parameters. Here, experimentalists may turn to lattice QCD to provide 
first-principles calculations of the factors required to coimect experiment with theory. 

In this proceeding, I will describe preUminary calculations of Bk on 2+1 flavor domain- wall 
fermion lattices. 



2. The Kaon Bag Parameter 

In the Standard Model, there are two sources of CP violation: direct and indirect. In indirect 
CP violation, decay into a disallowed number of pions is due a small amount of mixing between 
the and the K^. In order to use kaon data to constrain the unitarity triangle, we must disentangle 
the magnitude of kaon mixing, parametrized by the kaon bag parameter Bk . 

It is defined by 

= — ra^2 — ' (2.1) 

where Mk is the mass of the neutral kaon, fx is the decay constant of the kaon (given by its 
coupling to the axial current), and = {sd)L{sd)L is a four-quark operator coupling to left- 

handed quarks that changes strangeness by 2. 



3. Challenges 

In the chiral Umit, Bk contains only the operator of the form VV +AA: 

ff^vfAA = {srM^rnd) + (sYnrsd^srursd), o.i) 

which may be derived from L= {V —A) and parity symmetry, which forbids VA and AV terms. 
However, on the lattice we typically do not have the luxury of chiral symmetry; in this case, we 
must consider mixing with other chiral structures. 

Even worse, it may be shown in chiral perturbation theory that these operators are not propor- 
tional to mass in the chiral limit, while the desired operator is. If significant mixing occurs, we will 
see almost nothing of the continuum operator we desire. 

Domain- wall fermions offer a solution by allowing us to control the amount of chiral symmetry 
breaking by varying the size of the fifth dimension. The results of symmetry breaking are easily 
treated by the adding a residual mass. It is crucial that our residual mass be kept small, since the 
unwanted off-chirahty terms will contribute as 0{m^^^). 



4. Our Lattices 



Our action must be chosen to minimize chiral symmetry breaking while remaining computa- 
tionally fast enough to produce thousands of lattices. Since this means our lattices will be rather 
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Figure 1: Left: Pseudoscalar decay constant as a function of average valence quark mass (blue points mark 
physical/;;: and /k); Right: (PS|^'^'^"^|PS) (blue point marks physical kaon mass) 



coarse and our fifth dimension relatively small, we turn to improved actions to help lower the resid- 
ual mass. Adding a rectangle term to the action accomplishes this by smoothing the gauge field 
at short distances. We use the DBW2 rectangle coefficient (ci = —1.4069) and gauge coupling 
j3 = 0.72 in a gauge action of the form = — |- ((1 — 8ci ) Tr ?7piaq + ciTr ?7rect) ■ 

In summary, our lattices are 16^ x 32 x 8 and have bare sea quark masses = 0.04 and 
mi £ (0.01,0.02). We used the RHMC to generate 6000 trajectories and discarded 1000 for ther- 
malization. We select configurations separated by 50 for analysis, giving us a total of 100 sets of 
weak matrix elements for each light sea quark mass. Each set combines nondegenerate pairs from 
valence quark masses my G (0.01,0.02, ...,0.06). 



5. Preliminary Results 

The rho mass determines our lattice scale in the physical sea and valence quark mass limit ; it 
is consistent with the scale determined from the static quark potential. We find a^^ = 1.6(l)GeV. 

The residual mass is determined from the midpoint correlator, and agrees with the (negative) 
quark mass at which the pseudoscalar mass extrapolates to zero. We find mres = 0.0106(1), which 
is fairly large, being larger than our lowest light quark mass (0.01). 

The pseudoscalar mass then allows us to determine the light and strange quark masses. We 
find mi = 0.00171(21) and nis = 0.042(5). It is important that our strange quark mass be near our 
sea input (0.04 + mres), since we do not intend to extrapolate its value. 

The pseudoscalar decay constant enters in the normalization of Bk ■ It may be determined 
from a ratio of the wall-point pseudoscalar correlator to the wall-wall pseudoscalar correlator. 

_ 2(m, + mres) <;(0,0V2Mp 

where niq and mres are quark and residual mass, Mp is pseudoscalar meson mass, V is volume, and 
the '(^"s are correlators with superscripts denoting the source and sink operators (both Pseudoscalar 
in this case) and subscripts denoting source and sink shapes (point or wall). 
We fit fp to the chiral form 

fp = a + bm + cm\ogm^ (5.2) 
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Figure 2: Plateaus of B/> for highest (0.06) and lowest (0.01) average valence quark mass 

where m = + m^es- See Fig. |l| (left). 

The most naive way to derive the matrix element associated with Bk is simply to take the three- 
point correlator and divide by the wall-wall two-point correlator to remove the matrix elements 
associated with the wall sources: 



(PS|^|PS) =4Mj 



^wpw I'srci 



(5.3) 



'^ww (^src ; ^snk ) 

See Fig. (right). 

However, Bk may be derived at once without using the noisy wall-wall correlator by a clever 
use of wall-point correlators: 



Bf 



MjV 



'wpw 



2|K + /7Vs)2'^>r;(?src,0'^4^(^,^s„k) Zj 



727 



(5.4) 



where = 0.732 is the renormalization factor of the axial current computed on our lattices, 
ZffjT^ =0.93 is the renormalization factor of the three-point operator ^uT^^ ^ms = 1-02 
is the conversion factor from RI-MOM to MS-bar scheme. We take Zff/Z\ from previous calcula- 
tions on 2 flavor lattices and Z-^ from previous quenched calculations. 

We expect Bk to approach its asymptotic value far from the source and sink. Depicted in 
Fig. § are the heaviest and lightest valence quark masses {my = 0.01 and 0.06). The plateaus do 
not appear to have any unusual wiggles or trends. 

Bk is fit to a form including the chiral fits to fp and Mp: 



1/2 

{VS\ff\VS)=a + b—^\ogMl + cMl, 
Anfp 



(5.5) 



where Mp is a function of quark mass, but fp is taken in the chiral limit. See Fig. ||. 

Taking the valence=sea point (red on Fig. ^ for each set of lattices, we may extrapolate to 
the chiral limit for the light quarks (leaving the strange quarks at the physical strange quark mass). 
This yields a final value (See Fig. ^): 



Bif = 0.571(17) 



(5.6) 
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Figure 3: Pseudoscalar bag parameter as a function of average valence quark mass (blue point marks valence 
extrapolation to physical Bk, red point marks the valence=sea point) 




Figure 4: Two extrapolations of Bk to the physical point: red: extrapolates valence=sea to the physical 
point; blue: takes valence=physical, then extrapolates sea to the physical point 



Alternatively, we may take the valence=physical limit first (blue on Fig. taking advantage 
of the chiral fitting form, and then extrapolate linearly to the sea=physical limit. This alternate 
method yields a somewhat lower value of 

Bk = 0.556(24) (5.7) 



6. Conclusions 

Since neither extrapolation correctly accounts for nonlinear chiral terms, we take the average 
of the two and add a systematic error associated with the difference: 

B^ = 0.563(15)(21), (6.1) 

where the first error is systematic and the second statistical. In either case, the preliminary 2+1 
value of Bk is somewhat higher than the scale extrapolation from 2 flavor DWF values. See Fig. ^. 

In the future, we wish to improve this calculation with larger volumes, smaller residual masses 
and a full nonperturbative treatment of renormalization. Larger volumes will allow us to use longer 
plateaus to reduce statistical error and also diminish any finite-volume effects. Smaller residual 
masses will further diminish mixing with wrong-chirality operators. A correct treatment of the 
factors Zff and on 2+1 flavor lattices is necessary to complete this calculation, which will 
require NPR techniques. 
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Figure 5: Dynamical extrapolation of Bk to the physical point of 2+1 flavor lattices at 1.6 GeV (marked in 
bold black) compared to other recent values 
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